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Far-infrared Observations of the Very Low-Luminosity Embedded 
Source L1521F-IRS in the Taurus Star-Forming Region 
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ABSTRACT 

We investigate the environment of the very low-luminosity object L 152 IF IRS 
using data from the Taurus Spitzer Legacy Survey. The MIPS 160 fim image 
shows both extended emission from the Taurus cloud as well as emission from 
multiple cold cores over a 1° x 2° region. Analysis shows that the cloud dust 
temperature is 14.2 ± 0.4 K and the extinction ratio is Aiqq/Ak = 0.010 ± 0.001 
up to Ay ~ 4 mag. We find kiqq = 0.23 ±0.046 cm^ for the specific opacity of 
the gas-dust mixture. Therefore, for dust in the Taurus cloud we find the 160 /im 
opacity is significantly higher than that measured for the diffuse ISM, but not 
too different from dense cores, even at modest extinction values. Furthermore, 
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the 160 /iin image shows features that do not appear in the IRAS 100 /im image. 
We identify six regions as cold cores, i.e. colder than 14.2 K, all of which have 
counterparts in extinction maps or C^^O maps. Three of the six cores contain 
embedded YSOs, which demonstrates the cores are sites of current star formation. 

We compare the effects of L1521F IRS on its natal core and find there is 
no evidence for dust heating at 160 or 100 /xm by the embedded source. From 
the infrared luminosity Ltir = 0.024 Lq we find Lboi_int = 0.034 — 0.046 Lq, 
thus confirming the source's low-luminosity. Comparison of L1521F IRS with 
theoretical simulations for the very early phases of star formation appears to rule 
out the first core collapse phase. The evolutionary state appears similar to or 
younger than the class phase, and the estimated mass is likely to be substellar. 

Subject headings: circumstellar matter — stars: formation — stars: low-mass, 
brown dwarfs — stars: individual(L1521F IRS, L1527 IRS, IRAS 04368+2557 ) 



Introduction 



Investigation of starless cloud cores provides a way to study the initial conditions of star 
formation, where prestellar cores are those judged to be closest to gravitational cloud collapse. 



The L 1521F cloud core in the Taurus molecular cloud at 140 pc distance (iTorres et al 



(120071 ) and references therein) is a very dense (n ~ 10^ cm ^) and cold (T 9 K) cloud 



core containing 2 — 5M(n of gas detec t ed in NH.cj and other dense gas tracers (jCodella et al. 



19971 : lOnishi. Mizuno. fc Fukui Ill999l : ICrapsi et al. II2004I : IShinnaea et al.ll2004l ). 



Onishi. Mizuno. &: Fukui I (119991 ) singled out L1521F (their MC 27) as special in map- 
ping surveys of Taurus cores, and argued that the high density and infall asymmetry seen 
in the HCO^{'i — 2) line indicated a core in the earliest stages of gravitational collapse with 
a free-fall time scale of 10^ to 10'* yr. Consistent with the young age, IRAS observations 
indicated an upper limit of about 0.1 Lq for the luminosity of any possible embedded ob- 
ject. The presence of infall asymmetries i n other tracers and an advanced chemical clock 
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led other authors to similar conclusions (ICrapsi et al. I |2004J : iLee. Myers. &: Plumd 12004 



20071 ). Models of the cloud structure were consi stent with the picture of a dense 



cold core heated ext ernally by the interstellar radiation field ( lOnishi. Mizuno. fc Fukui 111999 



The Spitzer Space Telescope (IWerner et al. 1 120041 ) detected a very low-luminosity in- 
frar ed source in the L1521F prestellar core, a s ource too faint to ap pear in the IRAS cata- 
log (ITerebey et al.l l2005l : iBourke et al. 1 12006| ) . iBourke et al. I (120061 ) present submillimeter 
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and independent Spitzer data for the core. The source is detected at all IRAC and MIPS 
wavelengths ranging from 3.6 to 70 yum in wavelength. In this paper we further present 
160 fim data for the cloud. The source's low-luminosity suggests the object is substellar, 
and its similarity to class sources suggests it is extremely young. An object that pre- 
dates Spitzer, IRAM 04191-1-1522, is a l ow- luminosity ~ 0.1 L(7) c lass source in Taurus 



detected by ISO longward of 60 microns (lAndre. Motte. fc Bacmanrull999l: iDunham et al. 



2006 ) , that is now suggeste d as a very low- luminosity object (VeLLO) (Idi Francesco et al. 



2007: Dunham et al. 



2008|). However, IRAM 04191+1522 differs from L1521F by exhibiting 



both a molecular outflow and centimeter radio emission. Other compara ble objects detected 



by Spitzer include VeLLOs in the L1014 and L328 molecular clouds (lYoung et al.l 12004 



Lee et al. 


2005: 


Huard et al. 


2006) 



20061 ). The spectral energy distributi ons (SEP) are s i milar 



steeply rising with strong emission at 24 and 70 /xm in the infrared. IDunham et al. I (120081 ) 
present numerous VeLLO candidates from the Spitzer c2d survey. The low- luminosity of 
the VeLLOs suggests their masses are quite small. The presence of dense cores and infall 
asymmetries towards the better studied sources suggests some may represent a very early 
stage of star-formation, of accretion onto a substellar mass central object. 

Models of early cloud collapse define two distinct phases, the first protostellar core, 
when gas in the central region first becomes opaque at T ~ 100 K and heats up to T ~ 2000 
K at which point begins the second core phase, when molecular hydrogen dissociates and 
a hydrostatic object at stellar density forms. The properties of the first core depend sen- 



mass, 



it are roughly 0.01 — 


0.08 for 


(IBoss & Yorke 


1995 




Bate 


1998; 



Saigo. Tomisaka. fc Matsumoto 1120081 ). Numerical simulations suggest the model luminosity 
then shoots up to 1.0 Lq as the second protostellar core undergoes runaway collapse of a 
0.001 Mq, i.e. Jupiter mass object. The second core continues to accrete mass but other 
properties, such as luminosity, at times less than ~10,000 yr are still unknown, until the 
object becomes older and r nore m assive, and is detectable as a protostar or young brown 
dwarf. IStamatellos et al. I (120051 ) explore whether this phase would be observable, in the 
context of collapse models with high infall rates. Prior to collapse, the SED of the model 
cloud peaks near 190 fim due to heating by the interstellar radiation field. They find signif- 
icant additional heating occurs after the accretion luminosity exceeds 5.0 Lq, producing a 
peak wavelength as short as 90 fim during the class phase. In the context of t hese models. 



i t app ears VeLLOs are unlikely to modify the submillimeter SED of the cloud. lYoung et al. 



(120041 ) reach a similar conclusio n for the L1014 VeLLO s ource . Far-infrared wavelengths also 
problematic: previous work by iBoulanger fc Perault I (119881 ) using IRAS data found that 
the externally heated molecular cloud contributes a significant amount of 0.6Lq/Mq to the 
infrared luminosity. 
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To better understand the properties of L1521F IRS and its relation to early phases of 
star formation we analyze the MIPS 160 /im and IRAS 100 fim data near the source. The 
molecular cloud appears to dominate the emission at these wavelengths, and our goal is to 
better define the peak of the SED by determining the contribution of the embedded source 
to the flux. Using these data we also determine the properties of the L1521F cloud core at 
160 fim and relate it to other 160 /im cores in the region. 



2. Observations 



2.1. Spitzer Observations 



The Spitzer data near L1521F derive from the Taurus Spitzer Legacy Survey (jPadgett et al 



20091 ). a large 44 square degree map of the Taurus star- forming region at seven wavelengths 



ranging from 3.6 to 160 fim. The nominal spatial resolution of Spitzer is about 1.7" at 3.6, 
4.5, 5.8, and 8.0 fim with th e IRAC instrument, and is 6, 18, and 40" at 24, 70, and 160 fim 
with the MIPS instrument JPazio et al. 1120041 : iRieke et al. hooi l. Th e region near L1521F 
was observed on 2005-Feb-23 by IRAC and 2005-Mar-02 by MIPS. See lPadeett et all J2009h 
for an overview of the survey parameters and data processing. 

The 160 fim data obtained in fast scan mode mapping does not have enough redundancy 
to fill in completely all the gaps due to a dead readout and the intermediate gap between 
the array detector rows (MIPS Data Handbook V3.3 8/24/2007). Furthermore with only 
5 pointings per pixel, the effects of hard radiation hits and saturation translates into small 
regions without data. To deal with these gaps and to preserve as much as possible the diffuse 
emission, a 3 pixel by 3 pixel (15 " per pixel) median filter is applied to the image. The net 
effect is a slight redistribution of surface brightness (~ 15%) and smearing of the original 
beam from 40" to about 1' in size. To quantify the effect we measured the resulting spatial 
resolution in two ways. First we determined a FWHM of 52" for a theoretical 160 fim PRF 
that was processed using the same 3 pixel window filter. We note that most objects appear 
extended in the 160 fim mosaic. However we were able to measure an empirical average 
FWHM of 59" for five pointlike sources, which we take to be the effective spatial resolution 
of our 160 fim image. 

Table [1] gives the Spitzer position of the embedded infrared source, measured at 24 
fim to be least affected by extinction and hence most accurately detect the central source. 
Positions determined at shorter IRAC wavelengths are the same within the measurement 
uncertainties. The source is extended at IRAC wavelengths, hence the fluxes given in Table 
[2] depend on both the aperture and background determination. The flux density values 



are comparable to but differ from tliose in iBourke et al. I (120061 ) . Given tliat tlie source is 
extended we ascribe tlie difference as due to a difference in tlie pliotometric metfiod. 



2.2. IRAS Observations 



Images at 60 and 100 /im were obtained from tlie High- Resolution IRAS Galaxy Atlas 
(IGA), an atlas that includes the galactic plane, Orion, Taurus- Auriga, and p Oph regions 



(ICao et al.lll997l ). The coordinates of IGA images outside the galactic plane are given in 1950 
equatorial coordinates. In order to successfully precess and work with the images, several 
changes were made to the original FITS headers using the Goddard IDL astron library. First, 
a misspelled keyword RADECSYS was removed and replaced with the keyword RADESYS. 
Second, the keyword LONPOLE=0.0 was added. In addition, the reference pixels CRPIXn 
and reference positions CRVALn were changed to be at the image centers from their original 
off-image locations. This last crucial change appeared to minimize differences in various 
FITS software implementations of the CAR projection type. 

The IGA images were generated using the HIRES resolution enhancement algorithm and 
therefore have higher spatial resolution (2') tha n the native IRAS (3.8^ x 5.40, ISSA (4.30, or 



IRIS (4.3') survey images at 100 /zm wavelength ( ICao et al.lll997l : lMilville-Deschenes fc Lagache 



20051 ) . We constructed an empirical point spread function (PSF) using several sources in 
the Taurus survey region to find an effective spatial resolution at 100 /im of 140" x 110" 
(FWHM), about a factor of two times lower than that for the Spitzer MIPS 160 pm. data. 



3. Results 

3.1. Morphology of far-infrared emission 

In Figure 1, far-infrared images of dust emission near L1521F IRS show a striking 
similarity in the diffuse emission, indicating that the 160 and 100 /xm wavelengths are mostly 
tracing the same dust component in the Taurus cloud. The images are centered on the 
position of L1521F IRS, but throughout the 1.1° x 2.2° field the emission at 160 /im is 
extended with respect to the spatial resolution of 59" = 8300 AU = 0.04 pc at 140 pc 
distance. 

In the central 5' the L1521F cloud core is systematically brighter at 160 /im than at 



^http: / /irsa. ipac.caltech.edu/data/IGA/ 
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100 yum, a signature that the core is colder and denser than the surrounding cloud. In 
addition, the young stellar objects IRAS 04248+2612 and DG Tau appear as prominent 
point sources in the lower right corner and indicate regions of higher than average dust 
temperature. The ringing at 100 /im around the two point sources is an artifact of the 
HIRES processing. Other Taurus YSOs fall within the region shown (Fig. 1), but only the 
two named sources show significant emission at these wavelengths. 

Previous analyses of the diffuse emission seen by IRAS at 60 and 100 fim show the 
cirrus can be modeled as emission from dust associated with both the atomic and molec- 



ular components of the interstellar medium (ITerebev &: Fichlll986l : iBoulanger fc Perault 



19881 : IWood. Myers. &: Daugherty Ill994j : lAbergel et al.lll994l ). This extended component is 
variously referred to as galactic emission, infrared cirrus, or the diffuse interstellar medium. 
Boulanger et al. I (119881 ) model the diffuse IRAS cirrus as due to dust heated by the inter- 



stellar radiation field and find the dust is not in thermal equilibrium except at the longest 
(100 /im) wavelength. The morphology of the diffuse emission in the new MIPS 160 fim 
data shows it to be fundamentally the same as the IRAS 100 fim emission, albeit viewed at 
higher spatial resolution. 

However, emission at the shorter 60 fim wavelength necessarily implies the presence of 
warmer dust. This has long been understood as a consequence of the Planck function, where 
emission from cooler dust shifts to the Wien part of the spectral curve and becomes quite 
faint. Figure 2 shows extended emission along with numerous point sources in a 60 fim 
IRAS image; the 70 /xm MIPS image (not shown) shows the point sources at higher (~20") 
spatial resolution but has less sensitivity to the extended component. The morphology of 
the 60 fim diffuse emission in Figure 2 appears related to the 160 /im emission (Fig. 1) 
but not cospatial, as if representing emission from the surface layers of externally irradi- 
ated clouds. The embedded source L1521F IRS appears as a marginally detected point-like 
source in the center of the image. However, the object is found in the IRAS Point Source 
Rejects Catalog (i.e. detected b ut not reliable) as R04255+2644, which would explain why 



Onishi. Mizuno. &: Fukui I (119991 ) reported the source was not detected by IRAS. Apart from 
the faint point source, there is no extended emission in the central 5' of the IRAS image 
(confirmed by the 70 /xm MIPS image) comparable to that seen at 100 and 160 /xm (Fig. 
1) from the cloud core. Our focus in this work is the dense cold component and we do not 
further consider the 60 /im diffuse emission. 
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3.2. Analysis of extended dust emission 

Before determining the temperature and optical depth we first give expressions for the 
intensity. Assuming dust that is in thermal equilibrium at temperature T, then at frequency 
1/ the intensity is given by Jj, = (1 — e'^'')B^{T) where r,^ is the optical depth and B^{T) is the 
Planck function. In the limit of small optical depth, adequate for our purpose, the intensity 
becomes 

h^TMT). (1) 

From this expression we see that taking the ratio of the intensity at two frequencies leads 
directly to the dust temperature, if the frequency /wavelength dependent properties of the 
dust are known. In practice a range of dust temperatures will be present, but at a given 
frequency a limited range of temperatures contributes significantly to the intensity, as one 
'loses' cold dust whose emission falls on the Wien part of the blackbody curve. The contribu- 
tion of small stochastically heated grains can be neglected at these wavelengths; equilibrium 
thermal emission fr om big grains (> 0.25 /xm) dominates the emission at A > 70 — 80 /im 



f lLi fc Drainell200lh 



The optical depth can be expressed in terms of the specific opacity k^j in units of cra^g ^, 
the density of the gas-dust mixture p , and the path length L by, 

Tv = f^uPL 

= K^pninNH, 

= (TuNh, (2) 

and where p = 1.38 is assumed for the mean molecular weight per hydrogen atom, niH is 
the mass of a hydrogen atom, Nh is the hydrogen column density, and is the opacity in 
units of cm^ per H atom. The frequency dependence of the opacity is often expressed in 
terms of the power-law index /? where = i^oi^-'l^oY at some fiducial frequency z/q, and 
where (3 ranges from roughly 1.5 - 2 at infrared wavelengths. Based on their definitions, the 
relation between extinction and optical depth is given by A\ = 2.5log{e)Tx = 1.0857rA and 
moreover the extinction opacity + Kuahs includes both scattering and absorption. 

The scattering efficiency is specified by the albedo; at long wavelengths the albedo is zero, 
hence Ky = Huabs for no scattering, which is the regime that Eqn (1) describes. 



3.3. Dust temperature 



Figure 3 shows there is a strong linear correlation between the 160 and 100 pm images. 
This indicates the emission can be well modeled using a single dust temperature but varying 



- 8 - 



optical depth. The two regions immediately surrounding the sources IRAS 04248+2612 and 
DG Tau depart from the linear trend, which suggests there is internal heating by these 
YSOs. Table [3] gives the results of a linear fit to the data, after masking out the warm 
dust in a 6.5' x 6.5' square region near the two named YSOs. In addition, the 160 /xm 
image was smoothed to the approximate 100 /im resolution via convolution by the empirical 
IRAS 100 fim PSF, and the resulting images sampled at 7 pixels = 56", the approximate 
Nyquist frequency of the IRAS image. Figure 3 also shows evidence for a slowly varying 
background component, such as due to zodiacal or galactic emission, that appears as a 
nonzero intercept. As discussed in section 3.1, the origin of the extended emission has been 
extensively investigated for IRAS data. This background, which appears as a constant offset 
in our small field, is not obviously related to the Taurus cloud so we do not attempt to model 
it here. Our temperature determination, which is based on the slope, does not require us to 
make any assumptions about the background level. 

We find a dust temperature of 14.2 ± 0.4 K for this part of the Taurus cloud (region 
shown in Fig. 1), based on the slope of the fit (Fig. 3) and assuming (3 = 2 for the opacity 
power-law index. Color corrections were derived using the source emissivity from Eqn (1) 
and are modest (~ 10%) for both MIPS and IRAS data at this temperature. The value 
given for j3 = 1.5 indicates the uncertainty in the temperature if /3 di ffers from the nomina l 
value. An analysis using the completely independent COBE data (jhagache et al. 1998 ) 



shows excellent agreement; their Table 2 gives Td = 13.9 ± 0.2 K at a position in common 
with our Taurus image. 



3.4. Optical depth and opacity at 160 jixn 



To the extent that the dust temperature is constant, then the structure seen at 160 fim 
(Fig. 1) represents variations in optical depth, i.e. variations in column density. Using Eqn 
(1) and assuming a dust temperature of 14.2 K, then a constant factor 6.18 x 10^^ converts 
from intensity in MJy sr~^ to 160 fim optical depth (after background subtraction). Figure 
4 shows the resulting rieo optical depth map. Near the embedded YSOs, DG Tau and 
IRAS 04248+2612, the dust temperature is higher than assumed so the optical depth will be 
overestimated within 3' of the embedded sources. The optical depth can be underestimated 
in cold dense cores such as L1521F, a point that we return to. 

In order to convert optical depth to column density (or mass) requires the opacity; 

however the long wavelength opacity is uncerta in in molecular clouds, due to possible grain 

grow' t h, ice mantle format ion, and other factors ( Ossenkopf fc Henning IE994 : Stognienko. Henning &: Osser 



19951 : lEvans et al. Il200ll ). The situation is much better constrained for the diffuse ISM, 
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where iLagache et al. I ([19991) find a = 8.7 ± 0.9 10-^^(A/250 Mm) -^cm^ using COBE 
data, which appears compatible with IWeingartner &: Draine I (l200ll ) dust models. Scaling to 
the MIPS effective wavelength of 155.9 fim then the diffuse ISM opacity is 2.3 10~^^cm^ 
from COBE data, whic h is clo s e to the value of 2.0 10~^^cm^ for Ry = 3.1 dust from 



Weingartner fc Draine I (120011 ): iDraine I (120031 ). and references therein. In alternate units 



(using fj, = 1.38), the gas-dust opacity becomes 0.098 cm^ and 0.089 cm^ g~^, respec- 
tively. The opacity of dust-only is higher by a fact or of the gas to dust ra t io, th e value of 
which is Mgas/Mdust= 124 if using the models of IWeingartner fc Draine I (]200ll ); iDraine 
J2003h . 



To test whether the diffuse ISM opacity applies in Taurus we compare the 160 /xm 
extinction with extinction Ay, which gives an independent measure of the column density. 
There is a strong corres pondence (Fig. 5) between images of Ay at 4' resolution from 2MASS 
(IFroebrich et al. 1120071 ) and MIPS Aiqq, smoothed to the same resolution. The trend seen 
in Figure 6 is linear, although it exhibits more scatter than is seen in Figure 3 for the dust 
temperature. The value of the slope (Table Hj) is A^m/Ay = 0. 011 which is 2.6 times 
highe r than Amo/Ay = 0.00042 for the diffuse ISM (IDraine I I2OO3I : IWeingartner fc Draine 
20011). This comparis on is based on dust with Ry = 3.1, the kind found in the Taurus cloud 



(IKenyon et al. Ill994l ): however the predicted ratio hardly changes for Ry = 4.0 or 5.5 dust 



models. Comparing the observed and theoretical ratios we conclude the 160 /xm opacity in the 
Taurus molecular cloud is 2.6 higher than the diffuse ISM, namely a = 5.2 10~^^ cm^ and 
K, = 0.23 cm^ g^^ with estimated one sigma error of twenty percent. This indicates differences 
in Taurus cloud dust from diffuse ISM dust models, even in molec ular cloud regions with 
modest {Ay = 0.4 to 4 mag) extinction. Since the extinction maps of lFroebrich et al. I (120071 ) 
are in fact derived from near-infrared data we express our result relative to K band to find 
Aj^^iIAk = 0.010, using Ay/A^ = 9.29 from their assumed dust model. [Flaherty et al. 



(120071 ) report Ax/ Ax values for shorter Spitzer wavelengths. 



The results we present in Table H] are fundamentally based on an empirical comparison 
of the 160 /im intensity with Ay, hence we further assess the reli ability of the Ay determi- 



nation. In their analysis of a nearby region in the Taurus cloud lArce fc GoodmanI (119991 ) 
find significant differences between four different Ay methods. In that work they present Ay 
based on spectral types for 95 st ars; Figure 7 shows their Ay values versus 2MASS derived 
Ay from IFroebrich et al. I (120071 ). Comparison with a unit slope line shows good correspon- 
dence except for an offset, such that the 2MASS extinction is systematically 0.66 mag smaller 
than the Ay based on spectral type. The offset seems plausible given the 2MASS derived 
value represents an average over 4 ' resolution compared with the line-of-sight spectral value. 
We conclude that, apart from a nonzero offset, the Ay determination is approximately linear 
over the Ay = 0.4 to 4 mag range. 
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3.5. Changes in dust temperature 



As shown in Figure 3, the dust temperature near two young stellar objects, IRAS 
04248+2612 and DG Tau, differs from the average 14.2 K value. The luminosities o f 0.4 and 
6.4 Lq for IRAS 04248+2612 and DG Tau, respectively Jwhite fc Hillenbrand l[20oi ). further 
suggest that the dust "hot spots" are internally heated. We now consider the temperature 
near L1521F IRS, a source that is much lower in luminosity. Figure 3 shows that data 
points in the L1521F core follow the average trend, and hence show no signature of dust 
heating. Since the heating effect might be small, in Figure 8 we also plot the unsmoothed 
160 fxm versus the IRAS 100 fim data from Figure 1, sampled at the 160 /im resolution. 
The resulting increased structure at 160 /im leads to increased vertical scatter, as is seen. 
However notice that the points in the L1521F core all lie systematically above the other data 
points. This suggests a temperature effect leading to extra 160 fim emission from colder than 
average dust. We conclude there is no evidence for dust heating at 100 or 160 fim by the 
embedded infrared source in L1521F. We turn this statement into a flux density upper limit. 
Consider the peak emission value for L1521F, namely the highest 160 fim point in Figure 3 
for L1521F. If this data point were moved horizontally to the right, increasing its 100 micron 
intensity by about 3MJy sr~^ then we would plausibly conclude there is evidence for dust 
heating. We use our empirical 100 /im PSF to convert this intensity increment into an upper 
limit of 2 Jy for the 100 /zm flux density from circumstellar heated dust around L1521F IRS. 
The situation is different at shorter wavelengths where the emission is dominated by warm 
circumstellar dust rather than the cold cloud. The emission at 60 /xm shows this behavior 
(cf. Fig 2), where a near point- like source is detected, but not emission from the extended 
dense core. 

Previ ous studies at long er wavelengths demonstrate the presence of cold dust in L1521F. 
Moreover iKirk et al. I (120051 ) determine a characteristic value of 10 K from fitting the ob- 
served spectral energy distribution for a number of d ense cores. Temp eratures of some 



well determined Ta urus cores range from 8 to 12 K (IKirk et al. 1 120071 ). The models of 



Evans et al. I (120011 ) predict varying dust temperature such that the center of the prestellar 
core is the coldest region. The physical explanation is that the central region is partially 
shielded from the interstell ar radiation f ield b y the outer layers and thus can reach a lower 
temperature. Res ults fromlKirk et al. I (120071 ) give 9.3 ± 0.3 K for the dust in the L1521F 
core. In addition, ICodella et al.l (119971 ) find a gas temperature of T = 9.1 K in L1521F using 



ammonia observations. The gas and dust temperatures are similar and plausibly suggest 



that the gas and dust may be well coupled at the n 
cloud core. 



10^ cm ^ density found in the L1521F 
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3.6. Optical depth maps and 160 fim cold cores 

To accurately model the optical depth of dense cores requires knowing the dust tem- 
perature, which we know for L1521F but not for other cores since we lack longer wavelength 
data. Hence, we take a different approach to studying the cold core component. We a pri- 
ori assume that two different temperature components contribute to the 160 fim emission, 
namely T = 10 K cores and T = 14.2 K clouds. We furthermore use the IRAS 100 fim 
image to make our model of the cloud component. Thus we obtain an excess image that 
is sensitive to the col d cores by subt r actin g a suitably scaled IRAS 100 /xm image from the 



MIPS 160 /xm image. lAbergel et al.l (119941 ) used a similar method in the Taurus cloud, but 



they use 60 and 100 /im IRAS data, which is sensitive to somewhat warmer dust. 

Figure 9 shows the resulting core and cloud images, each converted to optical depth 
using the assigned temperature. Including 10 K dust means the optical depth is significantly 
higher (xl4) in core regions than in Figure 4. This implies that accurate optical depths 
for the cold cores depend on accurate dust temperatures. The dust temperature is known 
for L1521F, but has not been measured for the other cores, which adds significantly to the 
uncertainty. Despite the caveats, the core excess image is quite useful. In Figure 9 the core 
image consists of features that appear at 160 /im but not at 100 /im, while the cloud image 
looks like the 100 /im image. Due to the ringing artifact, we treated the regions around the 
two bright YSOs differently. We cut out regions around the two point sources and filled 
them with median values determined over a larger region. The effective smoothing this 
procedure introduces means the optical depth near the two bright point sources has higher 
uncertainty. Moreover, as for Figure 4, optical depths very close to the point sources (3') 
should be discounted, because those regions have warmer than average dust (see Fig. 3). 

One immediate use of the images is to identify dense cold cores, that show up more 
clearly than in Figure 4. The core image shows six distinct features that we identify as 
dense cold cores. Three out of six cores contain embedded YSOs, the previously described 



L1521F IRS, DG Tau, and 
extinction niap catalog by 



RAS 04248-1-26 1 2. Al l six cores have counterparts in the optical 



Dobashi et al. I (120051 ). The three cores containing YSOs were 



mapped by lOnishi et al. 1 Jl998h and designated as C^^O cores. Table [6] gives the 160 /xm 



cold core positions and cross identifications. 



3.7. Properties of the L1521F core 



To better understand our analysis we compare with previous studies of the L1521F core, 
starting with the extinction. The position of the infrared source is spatially coincident with 
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the peak of the extended emission at 160 /im; various studies show the core exhibits small 
scale structure so therefore the pea k extinction Ay should depend on spatial resolution. The 
extinction based on 2MASS data (IFroebrich et al. 1120071 ) has the advantage that it covers 
the entire Taurus cloud. However, the achieved spatial resolution of ~ 4' is relatively low 
and in high extinction regions the method begins to break down leading to values that are 
systematically too low. This is borne out for L1521F where Ay is the smallest we discuss: 
Ay is 3.7 toward the infrared source. The 1' resolution of our 160 /xm image is higher, and 
for our single dust temperature fit (Fig. 4) we find Ay = 8 at the peak position. The 
Ay drops to 6 if measured with respect to the image minimum, illustrating the uncertainty 
if a reasonable choice is made for a nonzero background level. We note that the assumed 
background level is not strongly dependent on image size, as the minimum in the L1521F 
image is close to the minimum found for the entire Taurus survey. 

However, our analysis shows L1521F is one of the regions that contains colder than 
average dust. The relevant result is from Figure 9 and the extinction increases substantially 
for our two temperature model (Fig. 9) to Ay = 43 mag, where the 10 K core contributes 
38 mag, and the 14.2 K cloud contributes 5 to the total extinction. With background 
subtraction the Ay drops to 41 ma g, a small ch a nge. A t yet higher (15") spatial resolution 



the trend of higher Ay continues: iKirk et al. I (120051 ) find a peak Ay = 130 mag based 



on 850 /im continuum data. In conclusion we find that analysis of the 160 //m emission 
provides a useful tool for determining Ay in high extinction regions, although very accurate 
dust temperatures are needed to fully utilize the method . 

Next we compute the mass of the L1521F core, and compare with previous studies. The 
mass M is given by M = pLA where A is the projected area of the cloud core. Recalling 
that = Tx, = then comparison with Eqn (2) shows the mass is related to the optical 
depth such that. 



M 



Atx/kx, 



0.61 Mo (0.23 cm^ g-'/Kx){D/UO pcfiVL^i^/S" x 8") J^ngo, 



(3) 



where ^pix is the pixel solid angle in square arcseconds and the sum is over pixels within 
the dense core, and after background subtraction. One source of error in the mass is the 
opacity, that we previously estimated to have 20% uncertainty. The opacity models of 
Ossenkopf &: Henning I (jl994j ) treat dust coagulation and ice mantle formation in dense core 
environments. In th eir study of dense core s 



0H2 and 0H5 from lOssenkopf &: Henning I (119941 ). where they favor 0H5 opacity (dust 



Evans et al. I (120011 ) consider dust opacities 



with thin ice mantles mn = 10^ cm ^ density gas). Both the 0H2 and 0H5 model opacities 
are consistent with our 160 /im opacity; specifically for Mg/Md = 124 the model 0H2 or 
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0H5 leads to 0.25 cm^ or 0.33 cm^ g~^, respectively, both falling within 2cr of our 
0.23 cm^ value. How to compare our value with submillimeter opacities is less clear, but 
assuming (3 = 2 and a gas to dust ratio of 100 instead of 124 gives Kgso = 0.01 cm? and 



'^1200 



0.005 cm? g~ 



values that are about a factor o f two lower than the 0H5 opacities 



but ar e essentially the same as the opacit i es ado pted by lAndre. Ward- Thompson, fc Motte 



( 119961 ) ; ICrapsi et al. I (l2004j ) ; iKirk et al. I (120051 ) . A second source of error is in finding the 
background level, that we estimate leads to 20% uncertainty in these extended objects. A 
third, and dominant, source of error is the temperature assumed for t he cold cores. W e 



adopted 10 K as a representative temperature for cold cores. However iKirk et al. I (120071 ) 
give T = 9.3 ± 0.3 K for the L1521F core; using this lower temperature leads to a factor 
of 1.9 increase in the optical depth and hence the mass. Presumably the 9.3 K figure is 
more accurate, but we also quote masses using the 10 K dust temperature to illustrate the 
uncertainty that may apply to other cold cores found within our field. 



The core mass determined by ICrapsi et al. I (120041 ) is 5.5 Mq from the 1.2 mm dust 
continuum; our measurement is 7.0, 3.7 Mq for 9.3 K, 10 K dust temperature, respectively. 
The mass is calculated using 136" radius aperture, minus a b ackground fr o m a sl ightly larger 
radius to mimic the on-the-fly mapping process. Similarly, iKirk et al. I (120051 ) find a core 
mass of 2.6 Mq from the 850 /xm dust continuum within 75" radius aperture; within the same 
aperture we find 2.8, 1.5 Mq for 9.3 K, 10 K dust temperature, respectively. Considering 
the difference in wavelength and techniques the agreement is good. However, in regions 
such as L1521F it is clear that an accurate dust temperature for the cold core is needed to 
determine an accurate mass. 

We point out that the same value of Kieo derived in modest Ay = 0.4 — 4 regions 
seems to work equally well for the L1521F dense core that contains both higher density and 
extinction. The opacity we determined for the Taurus cloud is higher than in the diffuse 
ISM, but we see no evidence for additional changes in the opacity within different molecular 
cloud environments. 



3.8. Infrared luminosity of L1521F 

Standard methods of determining bolometric luminosity are problematic for low-luminosity 
sources such as L1521F IRS because the emission longward of 100 fim is dominated by heat- 



i ng fr om the interstellar radiation field. To illustrate the difficulty, iBoulanger fc Perault 



( 119881 ) find an infrared luminosity per solar mass of gas of 1.8 Lq/Mq for the diffuse ISM 
and 0.6 Lq/Mq in molecular clouds; the YSO luminosity thus represents an excess above 
this 'background' luminosity. In the case of L1521F, prior to Spitzer, the source was success- 
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fully modeled purel y as a externally heated prestellar core (jOnishi. Mizuno. fc Fukui I Il999 



Crapsi et al. 1 120041 ) . more over whose emissi on longward of 100 /im could be fit by a sin- 



gle temperature greybody (IKirk et al. II2005I ). The usual definition of bolometric luminosity 
means to integrate luminosity over all wavelengths. Here we explicitly define Lhoi_int to mean 
just the portion of bolometric luminosity due to internal heating sources. To find Lhoi_int 
the ap proach taken by several authors has been to fit the source SEP based on a phy sical 
model (jPunham et al. II2006I : iBourke et al. II2006I ) and in this way lBourke et al. I (l2006l ) find 
a best fit value of 0.05 ± 0.02 Lq for L1521F IRS. Although the method is promising, there 
is a concern that the appli ed correction to the o bserved infrared luminosity is large and 
possibly model dependent. iDunham et al. I (120081 ) further develop and test this technique 
and suggest that Lhoi_int can be accurately estimated using the 70 /im luminosity. 

We adopt an alternative approach of reporting infrared luminosities based solely on 
the Spitzer data, so that sources can be compared with each other in a robust way. We 
then introduce SED templates to estimate the internal bolometric luminosity. One imme- 
diate problem is that numerical integration of steeply rising SEDs is sensitive to how the 
boundary is treated. To avoid this difficulty we adopt several extragalactic conventions 
(jSanders fc Mirabel Ill996l ) in part because galaxy SEDs share many similarities with em- 
bedded YSO SEDs. Specifically, Li,{X) = AnD'^fi, is defined to be the luminosity per unit 
frequency. The total infrared luminosity is then a weighted sum of 'ji'Li, over specified bands, 
where 7 is a coefficient of order unity. Several definitions of the total infrared luminosity. 



Ltir, have been proposed (see iHuynh et al. I (120071 ) and references therein). We adopt 
Ltir = 7iz^iv,,(8.0yum) + 72Z/L,,(24/im) + 73i/L,,(70yUm) + 74Z/Lj,(160/im) with 7 = [1, 1, 1, 1 ] 
for sim plicity, which is very similar to the 7 = [0.95,1.15,1,1] proposed by iDraine fc Li 
J2007h . 



One remaining task is to estimate the fiux density for L1521F at 160 /im since we 
detected only the cloud but not the embedded source at this wavelength. We do this by 
generating an SED template for class sources, using for our template the well-known class 
source L1527 IRS (IRAS 04368+2557) in Taurus. Because L1527 IRS is an extended 
source at 160 /zm we briefiy discuss its fiux density measurement. The fiux density is 54 
Jy using 64" radius aperture with background determined from a just-larger annulus. The 
choice of aperture and background is important, and to some extent arbitrary, which seems 
unavoidable. For L1527 IRS the measured HWHM size at 160 /im i s 50", meaning the 
64" aperture is slightly larger. To compare with KAO measurements of iLadd et al. I (jl99ll ) 
we use the background at 300"; our fiux density then increases to 76 Jy (es timated 20% 
uncert ainty), which compares well with the total fiux of 94 ± 33 Jy found by iLadd et al. 
( I1991I ). In our SED we use the 54 Jy fiux density as the value that better represents the 
circumstellar emission. Table [5] contains measured fiux density values for L1527 IRS while 
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Figure 10 shows the SED. 



The total infrared luminosity Ltir = 1.4 Lq is dominated by the 70 and 160 /xm fluxes, 
and moreover compares favorably with the IRAS result Lj^as = 1-6 Lq (IKenyon et al 



1990l ). Note that the different beam sizes, with IRAS being the larger (5'), likely account 



for an y differences. The source has been recently modeled by several groups. iFurlan et al. 
(120081) report L^ni ip .t = 1.8 Lq for the star plus disk luminosity from a model flt to the SED. 



Tobin et al. I (120081 ) present extensive modeling of the SED and IRAC images of this edge-on 
source. They determine Li,oi_int = 2.0 Lq from integrating an aperture-deflned model SED, 
and furthermore note that not all photons appear in the SED so the internal luminosity is 
in fact higher, leading to 2.6 Lq for the total value. In terms of the SED flt, the model 
undershoots the long wavelength part of the SED, which suggests that external heating by 
the ISRF is important, even for this relatively luminous YSO. To conclude, the Spitzer 
luminosity is Ltir = 1.4 Lq and based on the models, the internal bolometric luminosity 
ranges from ~ 1.9 Lq to 2.6 Lq which implies a multiplicative bolometric correction of 1.4 
to 1.9, respectively. 

Returning to L1521F IRS we flnd Lj-jr = 0.024 Lq (Tabled]) where we have estimated 
uLiyllGO fim) using the SED template. The estimated bolometric luminosity is then Li,oi_int = 
0.034 — 0.046 Lq, reflecting the uncertainty in the bolometr ic correction. This v alue is 
consistent with the lumin osity of 0.05 ± . 02 Lq found by the iBourke et al. I (120061 ) model 



and Lint = 0.03 Lq from iDunham et al. 
luminosity of the source. 



(120081 ). Our result conflrms the extremely low- 



3.9. Comparison of L1521F IRS with the class source L1527 



We continue our comparison of L1521F IRS with L1527 (IRAS 04368+2557), the one 
and only previously known class that falls within the boundary of the 44 square degree 
Taurus Spitzer Legacy Survey. Neither source was detected by 2MASS at 2.2 /im, consistent 
with their deeply embedded natures. However, the lower extinction at IRAC wavelengths (3.6 
- 8 /zm) and high sensitivity mean that Spitzer is sensitive to the outflow, showing reflection 
nebu losity with bipolar morphology toward both sources (IBourke et al. II2006I : iTobin et al. 
2008h . 



Figure 11 shows the two objects side by side and illustrates the difference in both spatial 
extent and source inclination. The symmetric morphology of L1527 IRS implies an edge-on 
z ~ 80 — 90° inclination, that is conflrmed by an 8 u m dip in the SEP (Fig. 10) that 



is due to deep silicate absorption (IFurlan et al. I l2008c iTobin et al. 1120081 ). By contrast 
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the morphology of L1521F IRS in both our data and that of iBourke et aL I (120061 ) finds 
the source near the apex of a conical, rather than a symmetric bi-conical structure at 3.6, 
4.5, and 5.8 /im; the fact that one outflow lobe is m uch brighter than the other in IRAC 



images implies the source has a modera te inclination (IWhitney et al. II2003I : iTerebey et al 



20061 ). The radiative transfer models of IWhitney et al. 1 J2003h and references therein show 
the importance of inclination on the SEP and images. B ased on our simulated images of a 
protostar with an opaque envelope (ITerebey et al. II2006I ) we estimate the inclination range 
is i ~ 50 — 70° for L1521F IRS. The smaller spatial extent as well as the moderate inclination 
of L1521F IRS further show that its luminosity is intrinsically low rather than low due to 
an inclination effect. 

The 1' scale bar in Figure 11 shows the approximate diameter of the infalling region, 
where the radius is given by r = at for a sound speed a = 0.2 km (appropriate fo r 
Taurus) and time t = 10^ yr (ITerebey. Shu, fc Cassen Ill984l : IShu. Adams, fc Lizano 1119871 ). 



The age of L1527 is not well known but recent results frorn the c2d Spitzer Legacy Survey 
give 1.3 X 10^ yr for the class lifetime (lEvans et al. II2008I ). The scale bar thus shows the 
approximate size of the collapsing cloud for class objects. If the VeLLO source were the 
same age, then in the context of low-mass star formation theory the mass of the collapsing 
region would be about 2Mq for the same size region (see section 3.7). 



4. Discussion 



4.1. Evolutionary status of L1521F IRS 



The interest of L1521F IRS and other VeLLO sources is where they fit into the evolu- 
tionary scheme of low-mass star-formation. Theory predicts an extended accretion phase, 
the embedded phase, during which a protostar grows in mass, and hence, luminosity, and 
displays the morphology of infalling cloud core (aka infall envelope), circumstellar disk, and 
outflow (ITerebey. Shu. &: Cassen 111984 IShu. Adams. &: Lizano 1119871 ). The embedded phase 
is represented by class and I sources, where class O's appear younger, based in part on being 
more de eply embedded but nonetheless having a ce ntral component capable of driving an 
outflow (lAndre. Ward- Thompson. &: Barsony Ill993l ). 



Although much lower in luminosity, VeLLOs share many similarities to known class 
sources, including a steeply rising SED and evidence for outflows. For the purpose of 
comparison we consider VeLLO sources separately from class sources, as their evolutionary 
status is still being deflned. For the VeLLO L1521F and class L1527, the cloud core 
masses are similar but the total infrared luminosity of L1521F IRS (0.024 Lq) is roughly 50 
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times smaller. The low-luminosity of L1521F IRS demonstrates it has a low central mass, 
regardless of whether the luminosity is dominated by L*, the central object luminosity, 
or Lace = fxGM^^M^:/ R^:, the accretion luminosity. Because the accretion luminosity is 
proportional to stellar mass, one interpretation for the low-luminosity of L1521F IRS is that 
it may be younger than L1527 IRS provided the mass infall rates are similar. Given the likely 
importance of initial conditions at this early phase it seems unlikely that the assumption 
of constant mass infall rate applies, s o it would be an overinte rpretation of theory to say 
that L1521F IRS is 50 times younger. IStamatellos et al. I (j2005[ ) predict luminosities above 
1 Lq an d would thus sugg e st Ve LLO sources are late class I phase. However, the infall 
rates in 



Stamatellos et al. 



(120051 ) are high compared with typical values for Taurus and 
the model duration is a relatively short 1.6 x 10^ yr to form an 0.5 Mq protostar. Using 
infall rates appropriate for Taurus we argue that the morphology and luminosity suggest 
that L1521F IRS could be younger or comparable in age to the class L1527, while the 
similar mass reservoir of the cloud cores suggests that over time L1521F IRS may eventually 
accrete enough mass to become a full fledged star. The cloud core has no well-defined edge, 
but contains at least 5 Mq, thus the source tests the idea that the final stellar mass is 



determined by the available mass reservoir and will approach a relatively constant fraction 



of the core mass 


(Motte. Andre. & Neri 


19981: 





nishi et al. 


19981: 


Testi & Sareent 


1998 


Johnstone et al. 


200 1|; 


Alves. Lombardi. & Lada 


20071: 


Enoch et al. 1 


2008|). 



At this point the statistics in Taurus suggest similar numbers of VeLLO and previ- 
ously known class sources in Taurus. To distinguish the two types we choose a luminosity 
boundary of Ltir = 0.5 Lq, which corresponds to roughly Lboi_int = 1-0 Lq. To motivate 
this choice we calculate that M^, = 0.1 Mq leads to an accretion luminosity Lace = 1-0 Lq, 
assuming a standard infall rate of 1 .6 x 10~^ M^yr"^, of which fraction = 2/3 falls onto a 
protostar of radius 1 Rq ( Fig. 7 of IStahler I (119881)') a nd is converted to accretion luminosity 
with fx = 0.62 efficiency (jShu. Shang. fc Leelll996l ). Further discussion of efficiencies ap- 



pears m 



Terebey et al. I (l2006l ). For sources whose luminosity is dominated by the accretion 
luminosity, the implied mass M^ = 0.1 Mq is close to the stellar-substellar boundary, so 
this choice of luminosity boundary distinguishes between objects that are stellar rather than 
substellar. 

In addition to L1521F and L1527, which lie within the Legacy Survey boundaries, there 
are two objects in the southern part of Taurus, namely the VeLLO IRAM 041914-1522 and 
the suggested transitional class /I source L1551NE (not to be confused with L1551 IRS5) 
dAndre. Motte. fc Bacmannlll999l : IPunham et al. 120061 : iMoriarty-Schieven. Butner. fc Wannier 
19951 ). If the star formation rate has been constant in Taurus then equal numbers implies 
equal time duration for the two types. Admittedly the numbers are so small that conclusions 
are premature, except to say that there is no evidence in Taurus that VeLLOs represent a 
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short-lived (10 x shorter) ph ase compared wi t h cur rent time scale estimates for class ob- 



jects. The survey results of iDunham et al. I (120081 ) present 15 new VeLLO candidates, so 



there is the promise of much better statistics in the near future. 

Existing theoretical simulations of the first protostellar core phase appear to be in- 



is relatively cold and produces litth 


i or no emission shortward of 70 um 


(Boss & Yorke 


1995; 


Bate 


1998: 


Stamatellos et al. 


2005: 


Saieo. Tomisaka, & Matsumoto 


2008|), whereas 



L1521F IRS is well detected at all IRAC wavelengths. The next phase, the formation of 
the second (hydrostatic) core, produces emission at the requisite IRAC wavelengths. How- 
ever, the existing simulations use high infall rates, and hence produce accretion luminosity 
L > 1 Lq, which is many times higher than the luminosity of L1521F IRS. A lower infall 
rate is plausible, hence the source L1521F IRS might represent the beginning of the second 
protostellar core phase. There is also the possibility that infall rates are variable. If the 
current infall rate is very low then the accretion luminosity can drop below for the cen- 
tral object. Theoretical evolutionary tracks are relevant for th i s case although still highly 







uncertain at young ages. The tracks presented in lBaraffe et al. I (l2002l ) suggest that 0.1 
is the maximum possible mass for L < 0.1 Lq if the age is less than 1 Myr. From this we 
conclude the current mass of L1521F IRS is likely to be substellar. 

Returning to Figure 11, we resume our comparison of the VeLLO and class sources. 
It is clear that L1527 and L1521F-IRS are different now (extent of nebulosity, luminosity, 
etc). However both sources could, in fact, have the same (substellar) mass and their dif- 
feren ces plausibly be explained by very different accretion rates onto the central protostar 



[e.g. 



Terebey et al. I (120061 )). If the accretion of L1527 were "turned off", or the accretion 
onto L1521F IRS "ramped up" then they could have the same SED, luminosity, CO outflow, 
and spatial extent. If the accretion rate were, for example, time variable then L1527 could 
represent a high-phase, and L1521F IRS a low-phase during class evolution. 

The data presented here do not offer a conclusive statement about the final mass of 
the protostar in L1521F IRS, whether it will remain a substellar object or continue to 
grow in mass, as the L1527 class source appears likel y to do. However, LI 52 IF IRS 



is embedded near the cent er of a collapsing cloud core (jOnishi. Mizuno. fc Fukui I Il999 



Lee. Myers, fc Plumd 120041 ). and it is likely to continue to accrete material (possibly at 
a variable rate) from an am ple mass reservoir until infall is terminated by one of possi- 



bly many ph y sical processes (jTerebev. Shu, fc Cassen 1 1 1984 IShu. Adams, fc Lizano 



1987 



Myers I l2008l : lEvans et al. I l2008l ). Given its location inside the cloud core, we suggest 
L1521F IRS will not remain substellar unless some condition, such as high angular mo- 
mentum or ejection, sucessfuUy terminates accretion soon. 
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5. Conclusions 

Using data from tlie Taurus Spitzer Legacy Survey, we analyze the 160 /xm far-infrared 
emission near L1521F to investigate the cloud environment and look for heating by the central 
object. The L 152 IF prestellar core contains a very low-luminosity embedded infrared source 
previously detected by the Spitzer Space Telescope. The source's low-luminosity suggests 
the object is substellar, and its similarity to class sources suggests it is extremely young. 

We investigate the influence of L1521F IRS on its natal cloud by studying the Spitzer 
MIPS 160 /xm and IRAS 100 /im data in a 1° x 2° region to look for evidence of dust heating 
at these wavelengths. In Figures 1, 3, and 5 we present the basic data and show (1) the 
emission is dominated by a uniform dust temperature of 14.2 K except near two bright YSOs 
in the field and (2) there is no evidence for any dust heating associated with L1521F IRS, 
suggesting that its impact on the dust continuum at A > 100 /im is minimal. Note however 
that at 70 /im the dust emission is dominated by the pointlike embedded source. 

The cloud dust temperature is 14.2 ± 0.4 K, based on the observed linear correlation 
between the MIPS 160 /im and IRAS 100 /im data. There is also a linear correlation with 
extinction up to Ay ~ 4, from which we derive Aiqq/Ak = 0.010 ± 0.001 for the extinction 
ratio and kiqq = 0.23 ± 0.046 cm^ for the specific opacity of the gas-dust mixture. 
The value of the opacity at 160 /im is consistent with dust models that incorporate dust 



coagu lation and/or ice mantle formation in dense core environments (lOssenkopf &: Henning 



19941 ). The 160 /im opacity is 2.6 times higher than that measured for the diffuse ISM, 
but is consistent with predictions for dense cores, which suggests that there are significant 
differences between dust in the Taurus cloud and the diffuse ISM even at modest extinction 
values. 

The MIPS 160 /im data show evidence for a cold core component that is not visible in 
the IRAS 100 /im data. The analysis produces an excess map at 160 /xm that is sensitive 
to colder regions, namely cold dense cores. In the subregion studied we find 6 cold cores, all 
of which have have counterparts in extinction maps or C^^O maps. By assuming 10 K dust 
temperature for the cold cores we compute the optical depth at 160 /im for both the cloud 
(14.2 K) and core (10 K) components. Three of the six cold cores contain embedded YSOs 
and thus indicate the cold cores are regions of current low-mass star formation. 

To measure source luminosity in an alternate way to model fitting, we introduce the 
Spitzer total infrared luminosity Ltir and make use of SED templates. For L1521F IRS, 
Ltir = 0.024 Lq and Lf,oi_int = 0.034 — 0.046 Lq, depending on the bolometric correc- 
tion. This result co nfirms the low source luminos ity of Lhni. m. = 0.05 L g^ determined by 



Bourke et al. I (l2006l ) and Lint = 0.03 Lq found by iDunham et al. I (]2008l ). In addition, we 
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estimate the source inclination is 50° — 70° based on the morphology of the extended emission 
in the IRAC images. 

To-date the Taurus Spitzer Survey contains the known class source L1527 IRS, as well 
as the VeLLO source L1521F IRS. The source L1521F IRS appears to be a substellar object 
that is at a similar or possibly younger evolutionary state than the class in L1527. The 
mass reservoir of roughly 2 — 5 Mq in the core suggests that L1521F is capable of eventually 
becoming a solar-mass type star. Comparison of L1521F IRS with theoretical models for the 
very early phases of star formation rule out the first core collapse phase. However it may be 
consistent with the early second collapse phase, defined as accretion onto a very low-mass 
hydrostatic object. Alternatively, if the infall rate is variable, and currently low, then the 
implied mass could be higher. However the low source luminosity is difficult to explain unless 
the source is substellar. 

This work is based in part on observations carried out by the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, Cahfornia Institute of Technology, 
under NASA contract 1407. S. Terebey warmly thanks the Spitzer Science Center, UCLA 
Dept of Astronomy, Caltech Dept. of Astronomy, and Harvard-Smithsonian Center for 
Astrophysics for their hospitality. M. Audard acknowledges support from a Swiss National 
Science Foundation Professorship (PP002-110504). 

Facilities: SSO, IPAC (IRAS). 
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Table 1. L1521F IRS Source Parameters. 



Name 


D 




RA 


Dec 




(pc) 


(Lo) 


(J2000) 


(J2000) 


L1521F IRS 1 


140 


0.02 


04:28:38.92 


+26:51:36.2 



Note. — Infrared luminosity from Table 5 and as defined 
in the text. Positions at other wavelengths are same as 
24 /im position, within measurement uncertainties. Position 
uncertainty at 3.6, 4.5, 5.8 /xm is 1.4", and at 8.0 /xm is 1.8". 



Table 2. Aperture corrected photometry of L1521F IRS. 



Telescope 


A 
(//m) 


PSF" 
(") 


Flux density** 
(mJy) 


a 

(mJy) 


Spitzer 


3.6 


1.66 


0.20 


0.04'^ 




4.5 


1.72 


0.35 


0.07'^ 




5.8 


1.88 


0.48 


0.10^= 




8.0 


1.98 


0.98 


0.2 " 




24 


5.9 


24. 


1.2 




70 


16. 


460. 


90. 


IRAS 


100 


118. 


<2000 





Note. — " Nominal PSF from Spitzer Observer's 
Manual (IRAC) or Data Handbook (MIPS). ^ IRAC 
photometry uses radius 2 pixel = 2.4" aperture and 2- 
6 pixel sky background. MIPS 24 uses PRF fitting. 
MIPS 70 uses radius 3 pixel = 12." aperture. '^Source is 
extended so flux depends on aperture and background. 
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Table 3. Taurus cloud dust temperature based on 160 and 100 fim data. 



Co 










^Icorr 


Tdust 














K 


7.94 


0.12 


0.259 


0.002 


2.0 


0.235 


14.2 










1.5 


0.235 


15.1 



Note. — Coefficients and uncertainty of linear 
fit to data 1100 = CO +C1*I160 shown in Figure 3. 

Corrected slope after applying MIPS and IRAS 
color corrections based on Eqn 1. 



Table 4. Comparison of 160 /im opacity with Ay. 





Aa 
^160o 








'^160 


MJy sr~^ mag~^ 


mag 








2 —1 

cm g 


18.1 


0.0018 


0.0011 


0.010 


5.2 10^2^ 


0.23 



Note. — Coefficients of linear ffi to data shown in Figure 5. " Estimated 
uncertainties are 10% although formal errors are much smaller. ^ Estimated 
20% uncertainties. Quantities at 160 /im are derived based on Ry = 3.1 dust 
model at wavelength V and with model parameters 4.896-E' — 22 = Cg^t 
extinction cross section {cm?/H), 1.870i? — 26 = M^ust per H nucleon 



(gram/H), 1.236^ + 02 = M„„jMd,,^t, 0.6735 = albedo flDraine 1 12003 : 



Weingartner &: Draine Il200ll ) 
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Table 5. Infrared luminosity of L1527 IRS and L1521F IRS. 



Source 




3.6 


4.5 


5.8 


8.0 


24 


70 


160 


TIR« 


L1527 IRS 


fuijnJy) 


1.2 


6.2 


11. 


8.2 


580. 


29000.^ 


54000.^ 






z/L.(A)(Lo) 


0.0006 


0.002 


0.003 


0.002 


0.043 


0.73 


0.59 


1.37 




uL^{X)/Ltir 








0.0001 


0.03 


0.53 


0.43 


1.0 


L1521F IRS 


fy{mJy) 


0.20 


0.35 


0.48 


0.98 


24. 


460. 








z/L,(A)(Lo) 


0.0001 


0.0001 


0.0001 


0.0002 


0.002 


0.011 


O.OIO'^ 


0.024*^ 




vLy{\)/LTiR 








0.009 


0.07 


0.49 


0.43^ 


1.0 



Note. — Fluxes from Table [2] or iPadgett et al.l (l2009l ). unless noted. More L1527 fluxes in 



Tobin et al. I (120081 ). Total infrared luminosity, see text. Source extended at 70 /im. Total flux 



density measured with standard 100" radius aperture, 120" - 140" radius background, 1.13 aperture 
correction. Estimated 20% flux density error. Source extended at 160 /xm includes diffuse cloud. 
Central flux density measured with 64" radius aperture, 64" - 96" radius background, background = 
170 MJy sr^^(i.e. extended cloud core), 1.38 aperture correction for T < 307^^ sources. Estimated 
20% flux density error. Estimated value based on L^/Ltir = 0.43, adopted from L1527. 
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Table 6. Cross Identification of Spitzer 160 /xm Cold Cores. 



ID 


R.A. 


Dec. 


/ 


h 


Ay Clump ^ 


C^^O Core2 




(deg) 


(deg) 


(deg) 


(deg) 


Name 


Number 


04:27:03+26:05:50 


66.76 


26.10 


171.84 


-15.68 


1211 P17 


21 


04:27:56+26:19:20 


66.98 


26.32 


171.81 


-15.38 


1211 P22 


22 


04:28:06+27:10:50 


67.03 


27.18 


171.17 


-14.78 


1211 P32 




04:28:39+26:51:40 


67.16 


26.86 


171.50 


-14.91 


1211 P16 


20 


04:29:13+26:14:20 


67.30 


26.24 


172.07 


-15.23 


1211 P24 




04:29:35+26:57:30 


67.40 


26.96 


171.57 


-14.69 


1211 P21 





Note. — The 160 /im peak position (J2000) using 15" pix els. 



^ Cloud name + clump name f rom 0.1° resolu t ion ca talog by lDobashi et al. 
^ Core name from Table [T] in lOnishi et al. 
resolution. 



(120051 ). 

fll998h and based on data with 0.05° 




04'^ 31 30 29 28 27 

Right Ascension 

Fig. 1. — Far-infrared images of dust emission in 1.1° x 2.2° field centered on L1521F IRS. 
The MIPS 160 iim (left panel) and IRAS 100 iim (right panel) images are strikingly similar, 
indicating the two wavelengths are mostly tracing the same dust component in the Taurus 
cloud. The L1521F cloud core (at image center) is brighter at 160 /im than at 100 //m, a 
signature that the core is colder and denser than the surrounding cloud. Two point sources 
are prominent in the lower right corner, especially at 100 fim, and indicate regions of higher 
than average dust temperature. The two bright sources are IRAS 04248+2612 and DG Tau 
from E to W. 
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04'^ 31"^ 30 29 28 27 

Right Ascension 

Fig. 2. — Extended dust emission is also seen across the IRAS 60 iim image but there is less 
correlation with the 160 /im image (Fig. 1). The 60 /im emission appears related but not 
cospatial, as if tracing warmer dust associated with the cloud surface. Field of view same as 
Figure 1. A faint but marginally detected source appears at the the position of L1521F IRS, 
exactly at the image center. Most point sources seen at 60 /im do not appear at longer 
wavelengths. 
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Fig. 3. — There is a strong linear correlation between the resolution matched MIPS 160 //m 
and IRAS 100 iim emission. Solid hne shows a linear fit to the data while symbols show 
regions surrounding specific objects. Excess 100 /im emission indicates warm dust around 
the YSOs IRAS04248+2612 and DG Tau. No evidence of warm dust is seen towards L1521F 
IRS, suggesting the source luminosity is too low to provide significant dust heating. 
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04^ 31 30 29 28 27 

Right Ascension 



Fig. 4. — Optical depth map at 160 /iin using constant dust temperature (using Table [3]). 
Solid lines show 0.0032, 0.0044, 0.008 contour levels. If measured with respect to image 
minimum (0.002 at 04:29:12.7 25:49:13 position) then the contours correspond to 6Av = 
1,2,4, respectively (using Table Hj). Optical depths are less accurate at the highest contour 
level because the constant temperature assumption breaks down in these regions (see text). 
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04"^ 31 30 29 28 27 

Right Ascension 



Fig. 5. — Compariso n of extinction de t ermin ed from 2MASS with that from 160 fim data. 
Left panel is Ay from lFroebrich et aL I (120071 ): right panel is extinction Aiqq (see Fig. 4) but 
smoothed to have matching 4' resolution. 
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Fig. 6. — There is a linear correlation between Ay and Aiqq up to Ay ~ 4 for the spatial 
resolution matched data. 
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1 2 3 4 5 

Av_2moss(nnog) 

Fig. 7. — Comparison of two different Ay metfiods shows good agreement with unit slope up 
to Av ~ 4. However the two methods have a systematic 0.66 mag offset, and above Ay ~ 2 
the scatter increases. 
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20 



30 



40 



(lOOAim) MJy/sr 



Fig. 8. — Using the unsmoothcd 160 /im data more clearly shows the 160 fim excess from the 
cold dust in the L1521F core. Regions near the two embedded sources have been excluded. 
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04"^ 30 29 28 27 



Right Ascension 

Fig. 9. — Optical depth greyscale maps at 160 /xm for core component (left) and cloud 
component (right). The 160 iim excess map (left) traces regions that differ from the average 
behavior (right) . The emission has been converted to optical depth for cores (left) assuming 
Td = 10 K and for the cloud (right) using = 14.2 K. Adding the maps gives the total 
extinction. Optical depth contours start at 0.03 with 0.03 increment. 
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Fig. 10. — Comparison of class and VeLLO spectral energy distributions. The top SED 
shows the central point source of L1527 IRS from 3.6 fim to 160 fim. The dip at 8.0 fim 
i s associated with 10 (j,m silicate absor ption that is due to the edge-on source inclination 
JPurlan et al. Iboosi : iTobin et al.lboosi ). The overall SED of L1521F IRS is lower indicating 
a lower luminosity. Both sources have steeply rising SEDs showing they are deeply embedded. 
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Fig. 11. — Comparison of class and VeLLO images at IRAC 4.5 fim wavelength shows 
scattered light nebulosity toward both. Scale in green is 1' = 8400 AU. North is up, east to the 
left. (Left) Image of L1527 IRS edge-on system shows a disk oriented N-S and outflow cavity 
E-W in scattered light. (Right) L1521F IRS lies at apex of conical nebulosity, suggesting 
an outflow seen at moderate inclination. The smaller extent of the VeLLO nebulosity is 
consistent with its lower luminosity. 



